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Resonant and inelastic Andreev tunneling 
observed on a carbon nanotube quantum dot 

J. Gramich, A. Baumgartner,* and C. Schdnenberger 

Department of Physies, University of Basel, Klingelhergstrasse 82, CH-4056 Basel, Switzerland 

We report the observation of two fundamental sub-gap transport processes through a quantum dot 
(QD) with a superconducting contact. The device consists of a carbon nanotube contacted by a Nb 
superconducting and a normal metal contact. First, we find a single resonance with position, shape 
and amplitude consistent with the theoretically predicted resonant Andreev tunneling (AT) through 
a single QD level. Second, we observe a series of discrete replicas of resonant AT at a separation of 
~ 145 /xeV, with a gate, bias and temperature dependence characteristic for boson-assisted, inelastic 
AT, in which energy is exchanged between a bosonic bath and the electrons. The magnetic field 
dependence of the replica’s amplitudes and energies suggest that two different bosons couple to the 
tunnel process. 

PACS numbers: 74.45.-he, 73.23.-b 


Electron transport in nanostructures with supercon¬ 
ducting (S) and normal metal (N) contacts have attracted 
considerable attention recently, because of the discovery 
of new transport processes and exotic quantum mechan¬ 
ical states of matter. While S contacts have been used in 
bias spectroscopy for many years, transport in nanostruc¬ 
tures at energies below the superconductor’s energy gap 
became relevant only recently, for example in the search 
for solid-state versions of Majorana Fermions [1], Cooper 
pair splitting [2-5] , or the study of Andreev bound states 
(ABSs) [6-9]. 

The most fundamental low-energy transport process 
between an N and an S reservoir is Andreev reflection, 
in which an electron from N can only enter S by form¬ 
ing a Cooper pair with a second electron of opposite spin 
and momentum (for a standard s-wave superconductor). 
This process is slightly more complicated if the two reser¬ 
voirs are connected by a quantum dot (QD), in which the 
charging energy suppresses a double occupation. The re¬ 
sulting resonant Andreev tunneling (resonant AT) is il¬ 
lustrated in Fig. la as an intuitive sequential tunneling 
process. At zero bias the two electrons forming a Cooper 
pair tunnel through the QD at the same energy, which 
leads to a peak in the differential conductance G. At a 
finite bias, resonant AT is allowed only for electrons with 
energies aligned to the electrochemical potential fis of S, 
i.e., only if the QD resonance (/igo) is aligned to //s? see 
Fig. la. Resonant AT can be identified by the distinc¬ 
tive resonance line shape with a sharper decay than when 
tunneling into a normal reservoir [10-12]. 

An additional transport process is boson-assisted tun¬ 
neling, in which bosons from the environment are ab¬ 
sorbed or emitted in an inelastic tunneling process. This 
mechanism leads to replicas of the elastic tunnel process 
at higher energies. For example, in carbon nanotubes 
(CNTs) with normal metal contacts, phonons can result 
in additional resonances [13] and even in the suppres¬ 
sion of the current for strong electron-phonon couplings 
(Franck-Condon blockade) [14]. Also the electromagnetic 


environment can provide energy for photon-assisted pro¬ 
cesses, resulting in discrete QD resonance side bands 
[15, 16]. Phonon [17] and photon [18] induced discrete 
replicas of elastic AT were predicted recently for QDs 
with an S-contact. Such inelastic AT is illustrated in 
Fig. lb: a Cooper pair can only be formed by two elec¬ 
trons of opposite energy with respect to the chemical 
potential of S, /is- If the QD level is aligned at a positive 
energy, both electrons traverse the QD at the same en- 
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FIG. 1: (Color online) (a) Schematic of resonant AT. Two 
electrons tunnel through the QD into S for /xqd = /is and 
/iN > /iQD. (b) Schematic of inelastic AT. Two electrons tun¬ 
nel through the QD into S at /igo = /is + f for /iN > /igo 
by emitting multiples of the energy hwh to the environment, 
(c) SEM image of a representative device and schematic of 
the measurement setup, (d) Differential conductance of QDl 
as function of Fsd and Fbg- (e) Detailed conductance map 
of the region indicated in (d). The resonant AT line and the 
inelastic AT are labeled R and IE, respectively. The dotted 
lines point out the Coulomb blockade diamond edges and ABS 
an Andreev bound state. 
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ergy, and, to form a Cooper pair, one relaxes the energy 
fkJb to the environment with a bosonic excitation spec¬ 
trum. As is evident from Fig. Ic this condition is met 
at /iQD — fJ'S = ^ N), in contrast to twice this 

value in boson-assisted processes in devices with normal 
contacts. Neither resonant nor inelastic AT have been 
observed in QDs before, probably because both require 
the QD resonance width to be much smaller than the 
superconducting energy gap, A. 

Here we use a CNT quantum dot contacted by a su¬ 
perconducting Nb and a normal metal contact to inves¬ 
tigate a discrete subgap resonance fully consistent with 
resonant AT. In addition, the resonance has multiple, 
very sharp replicas due to inelastic boson-assisted AT. 
Nb shows a large energy gap in CNT devices [9, 19], in 
our devices up to 1 meV, which allows us to perform bias 
spectroscopy in the regime of the QD life time broaden¬ 
ing being considerably smaller than A, which is crucial 
to identify these processes. 

Figure Ic shows a false color scanning electron mi¬ 
croscopy (SEM) image, including a schematic of the mea¬ 
surement setup. The CNTs were grown by chemical va¬ 
por deposition (CVD) on a Si++/Si02 substrate used as 
a backgate. Subsequently the surface is treated by an 
rf-induced hydrogen plasma, similar to [9, 20]. In this 
process about 30% of the CNTs are etched away, which 
probably serves as a selection process for defect-free 
clean CNTs. On a suitable CNT we fabricate a central 
^ 200 nm wide and ^ 2 mm long Ti/Nb (3nm/40nm) 
contact (S) and two normal metal Ti/Au (5 /65 nm) con¬ 
tacts (N) on either side of S at a distance of ^ 300 nm, in 
a Cooper pair splitter geometry [5] with sidegates (SGs) 
for individual electrical tuning of the two CNT sides. 
Between the contacts two separate QDs form, but no 
signals could be found that depend on both QDs [21]. 
Here we focus solely on experiments on QDl. The ex¬ 
periments were carried out in a dilution refrigerator at 
a base-temperature of ^ 110 mK. In Fig. Id the differ¬ 
ential conductance G through QDl is plotted as a func¬ 
tion of the backgate voltage, Vbg^ and the bias applied 
to S, Esd- We find very clear Coulomb blockade (CB) 
diamonds with an addition energy of ^ 6.6 meV, a level 
spacing of ^ 1.2 meV, and a well-defined superconducting 
gap of ^ 1.2 meV, which separates the CB diamonds in 
a characteristic way due to the gap in the single-particle 
density of states in the S contact [21]. This gap is reduced 
monotonically with a magnetic field B applied perpen¬ 
dicular to the CNT and the substrate and vanishes at 
H « 1.5T [21]. Fig. le shows G plotted in the region 
indicated by the rectangle in Fig. Id. Here we find up to 
7 weak but very sharp parallel resonance lines with the 
same positive slope as the CB diamonds. Similar fea¬ 
tures appear in all other CB diamonds [21], independent 
of the QD charge state. The average spacing between the 
lines is ^ 145 ± 30 //eV, with some resonances showing 
significant deviations from this value. Only one of these 
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FIG. 2: (Color online) Resonance line shapes: G as function 
of Vbg (points) at zero bias (a) in the normal state at H = 5 T 
and (b) at zero field in the superconducting state. The blue, 
red and dotted green lines are best fits to the data using 
the expressions for a Breit-Wigner (BW), resonant AT and a 
thermally broadened QD resonance, respectively. 


parallel lines (R) crosses the entire transport gap, while 
the other lines (IE) end at a finite sub-gap bias. Pro¬ 
nounced negative differential conductance (NDC) values 
occur between the conductance maxima. An additional 
structure arises due to an Andreev bound state (labeled 
ABS in Eig. le), which leads to some amplitude mod¬ 
ulation, but will not be discussed in more detail here. 
Line R we identify with resonant AT. To demonstrate 
this, we compare the line shape of the zero-bias CB reso¬ 
nance measured in the normal state at H = 5 T shown in 
Eig. 2a to the resonance at zero field, i.e. in the supercon¬ 
ducting phase of S, which is plotted in Eig. 2b. Eor the 
following fits we added an identical background deter¬ 
mined by the data points far off resonance (not shown). 
In the normal state we expect a Breit-Wigner (BW) line 
shape due to life-time broadening, which is described by 
G(AVbg) = X ’ /4 , with F = ri+r2 and the QD 

level detuning AE = —ea{AVBG — ^g)’ where a is the 
lever arm of the backgate to the QD and the posi¬ 
tion of the resonance. This expression fits the observed 
line shape in Fig. 2a very well (blue curve) for the tunnel 
coupling parameters Fi ^ 9.0/reV and r 2 ~ 96.5/ieV. 
As a comparison, we also plot the best fit using the 
expression for thermally broadened resonances (dotted 

green line), G = x 4 /cbT cosh~^ which does 

not describe the data well. For resonant AT in the limit 
of non-interacting electrons the expected line shape is 

G = ^ which also deviates 
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FIG. 3: (Color online) Temperature dependence of resonant 
and inelastic AT. Differential conductance G as function of 
Vbg at (a) VsB = OmV and (b) Vsb = —0.5 mV for the 
three indicated temperatures. The black arrow points out the 
resonant AT line, the red arrows the boson-absorbing and the 
blue arrows the boson-emitting resonances [21]. The inset in 
Fig. (a) shows the peak amplitudes of the features indicated 
in the main figure. 


strongly from the data. In contrast, the resonance at 
zero magnetic field plotted in Fig. 2b is described best 
by the expression for resonant AT: the measured con¬ 
ductance values decay faster away from the maximum 
than in the normal state. The extracted values for the 
tunnel coupling using the BW expression are a factor of 
3 — 7 smaller than in the normal state, while the values 
extracted from the resonant AT expression, Fi « 8.4 /ieV 
and F 2 ~ 66.5 yueV, are very similar to the ones in the 
normal state. We find similar values also in gate sweeps 
at a small bias [21]. From the line shapes and the obser¬ 
vation that this resonance occurs when the QD level is 
aligned to the electrochemical potential of S, we conclude 
that resonance R running through the full energy gap of 
the superconductor is due to resonant AT. 

The replicas (IE) parallel to the resonant AT line (R) 
found in Fig. le we attribute to the emission and ab¬ 
sorption of energy from or into a bosonic reservoir in 
an inelastic Andreev tunneling process. The bias condi¬ 
tion for this process (see Fig. lb) is /is < /^qd < Mn 
[ 17, 21], which results in the lines ending at a finite sub¬ 
gap bias, as observed in Fig. le. Following the pre¬ 
dictions of Ref. [17] for phonons. Fig. 3 shows a series 
of backgate sweeps obtained at different temperatures. 
The zero-bias curves in Fig. 3a are symmetric with re¬ 
spect to the central resonant AT peak indicated by the 
black arrow. The amplitude of this peak is reduced con¬ 


siderably with increasing temperature, while the width 
increases because of side peaks (red arrows) emerging as 
weak shoulders between 0.5 K and 1 K. The amplitudes of 
the indicated features are plotted in the inset of Fig. 3a 
for the three temperatures, which shows that resonant 
AT is reduced with increasing temperature (^ 1 //cbT) 
due to the thermal broadening of the Fermi functions in 
the normal metal contact. In contrast, the absorption 
side peaks increase in amplitude at higher temperatures 
due to the thermal population of the boson states [28], 
in good qualitative agreement with the predictions for 
phonon-assisted AT [17, 29]. 

Backgate sweeps at the bias Vsd = —0.5 mV are plot¬ 
ted in Fig. 3b for the same temperatures, which again 
shows the reduction of the resonant AT amplitude (black 
arrow) and the onset of boson-assisted inelastic AT at 
more negative gate voltages (red arrows). At more pos¬ 
itive voltages, i.e. /is < Mqd, resonances occur with 
slightly decreasing amplitudes and essentially constant 
widths (blue arrows) [21]. In addition, we observe pro¬ 
nounced negative differential conductance (NDC) be¬ 
tween the peaks. All these findings are in very good 
qualitative agreement with inelastic AT in which the ex¬ 
cess energy of the tunneling electrons is emitted into a 
bosonic bath [17]. In particular, the resonance width 
in this process is mainly determined by the tunnel cou¬ 
plings, which we do not expect to change greatly with 
temperature. NDC is also expected, since the resonance 
condition is determined by the QD level position and not 
by the Fermi energy of the normal metal lead, which leads 
to peaks in the current and a peak-dip structure in the 
differential conductance. 

We now investigate the resonant and inelastic AT as 
a function of an external magnetic field B applied per¬ 
pendicular to the substrate plane. Figure 4a shows G 
as a function of Vsd and B at the backgate voltage at 
which the elastic AT resonance (R) crosses Vsd = 0. 
We find that the energy gap detected by the QD shrinks 
monotonously up to ^ 0.3 T and is then roughly con¬ 
stant up to ^ 1 T and diappears around 1.5 T [21]. Here 
we focus only on the sub-gap features: the resonant AT 
line R is essentially unaffected at fields below B « 0.7 T 
and splits at higher fields due to Zeeman shifts of the 
resonances. These do not affect the following analysis 
since we extract all values directly from Coulomb dia¬ 
mond experiments [21]. We can visualize the transition 
from resonant AT to a Breit-Wigner characteristics with 
increasing field by plotting the extracted tunnel param¬ 
eter, F = Fi+F 2 , asa function of H, which is shown 
in Fig. 4b. We assume that the tunnel couplings stay 
roughly constant with increasing B [21]. While up to 
B ^ IT the expression for resonant AT reproduces nicely 
the large-field values of the BW fit (light red band), it 
overestimates the coupling by more than a factor of 2 
at high fields. At intermediate fields between 1 T and 
3T both line shapes deviate considerably from the ex- 
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FIG. 4: (Color online) Magnetic field dependence of resonant 
and inelastic AT. (a) G as a function of the bias Vsd and the 
external magnetic field B at a fixed backgate voltage where 
the resonance (R) crosses Rsd = 0. (b) Tunnel coupling T = 
Fi +r 2 of resonance (R) extracted from the resonant AT and 
the Breit-Wigner expressions as a function of B. The fits are 
obtained from colorscale images always at the same relative 
position [21] (c) Energy and (d) amplitude of the inelastic AT 
peaks (IE) as function of B. The resonances are labeled in 
Eig. c for increasing energy. 


pected values, which suggests that normal electron tun¬ 
neling and AT co-exist at fields where the visible trans¬ 
port gap is reduced to zero. 

Figure 4c and 4d show the magnetic field dependence 
of the position (energy) and the peak amplitude of the 
inelastic AT lines, respectively, extracted from CB spec¬ 
troscopy [21]. We plot the position of the conductance 
maxima relative to the resonance R at negative and pos¬ 
itive bias up to IT and label the lowest resonances in 
energy in Fig. 4c. While the energy of resonance lEl is 
essentially constant in this field interval, the positions of 
IE2, IE3 and IE4 are all reduced and scale similarly with 
increasing field, but not linear with the energy gap [21]. 
For the latter resonances the spacing stays roughly the 
constant. Even more pronounced is the difference in the 
field dependence of the resonance amplitudes, plotted in 
Fig. 4d: the amplitude of lEl is independent of B within 
experimental error, but the amplitudes of IE2 and IE3 
decay continuously on the scale of 1 T. 

Inelastic AT is mediated by the absorption or emission 
of bosons. In our system, three types of bosons might be 
responsible for these sub-gap processes: (I) phonons, e.g. 
mechanical oscillations in the CNT [17], where the longi¬ 
tudinal acoustic modes in the doubly-clamped CNT, or 
the CNT squash mode [30] can yield quantized energies of 
the observed energy scale. (II) plasma modes in the mil¬ 
limeter scale superconducting contact [31]. (Ill) photons. 


i.e. electromagnetic modes of the resonator formed by the 
inductance and capacitance of the S-contact, which are 
damped by ohmic dissipation in normal metals. Phonons 
naturally account for the temperature dependence of the 
inelastic AT lines, but it is not straight-forward to ex¬ 
plain a magnetic field dependence. Photons would, at 
least qualitatively, account for the field and temperature 
dependence, since the kinetic inductance Lk of S diverges 
with decreasing energy gap, which results in a reduced 
resonance frequency oc (LkC)”^/^ {C is the capacitance 
to the backgate) [32] . The fact that we find qualitatively 
different magnetic field characteristics for the lowest en¬ 
ergy inelastic AT peak compared to the resonances at 
higher energies might indicate that two different bosonic 
baths are coupled to our QD. 

In summary, the large superconducting transport gap 
found in our Nb-contacted CNT QDs and the sharp QD 
resonances allow us to identify resonant (elastic) and in¬ 
elastic Andreev tunneling in a QD-superconductor struc¬ 
ture. In addition, we also demonstrate the impact of the 
energy gap on the CB diamond structure. The tempera¬ 
ture dependence of the inelastic replicas of resonant AT is 
consistent with bosonic excitations that open additional 
transport channels. However, from our experiments the 
nature of the bosons is difficult to assess. The magnetic 
field dependence might even hint at the possibility of two 
different bosonic systems coupling to the QD. Our experi¬ 
ments demonstrate that, in contrast to normal metal sys¬ 
tems, such excitations can become the dominant trans¬ 
port mechanisms in S-QD systems, while a smearing of 
the discrete resonances is a possible origin of the heav¬ 
ily discussed ’soft gaps’ [33] in large-gap superconduc¬ 
tor nanostructures. In addition, the coupling to bosonic 
reservoirs is also expected to result in replicas of other 
sub-gap features, for example ABSs [34]. Similar devices 
with engineered bosonic environments could shed more 
light on the nature of the observed processes and lead 
to hybrid quantum systems with a controlled coupling 
to the superconductor. For example, novel techniques 
allow the fabrication of suspended CNTs with supercon¬ 
ducting and other contact materials [35], which should 
give access to well-controlled CNT phonons [36]. Using 
the same technique, a radio-frequency cavity can be cou¬ 
pled to the CNT QD, which results in a well-controlled 
electromagnetic environment with discrete modes [37]. 
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SUPPLEMENTARY MATERIAL 
QUASI-PARTICLE TUNNELING IN AN N-QD-S SYSTEM 

We briefly discuss quasi-particle tunneling in a S-QD-N structure. The corresponding experimental results can 
be found in Fig. Id and le of the main text and in Fig. S8. Figure S5 sketches the differential conductance for a 
quantum dot (QD) connected to a normal metal (N) and a superconducting contact (S) with the bias applied to S (see 
Fig. la main text). Figure S5a shows the relative positions of the electrochemical potentials /iN, = —eVsD + Mn 
and fiQB for selected positions of the expected charge stability diagram. The tips of the Coulomb blockade (CB) 
diamonds mark the onset of quasiparticle tunneling, separated in bias by 2A/e. The tips are shifted in gate voltage 
with respect to each other due to the capactive coupling of the leads to the QD by AV^ = 2A/(e/3_), with /3_ = 
—Cg/{Ctot — Cs) the (negative) slope of the /iN = /^qd CB resonance. Cs, Cn, Cg, Ctot = Cs + Cn + Cg are the S, 
N, gate and total capacitances, respectively. If one considers only single electron transport through the system and 
neglects superconducting correlations and charge dynamics on the QD, the qualitative charge stability diagram can 
be understood by a simple resonant tunneling picture. The current is given by[l] 

/ (X / dE^^{E)&s{E + eVsD) ■ Tqb{E) ■ [ME) - fs{E + eVsD)]. (1) 

The resulting differential conductance G = d//dUsD is plotted in Fig. S5b which reproduces the qualitative arguments 
made in Fig. S5a. Here we chose the density of states (DOS) ^^{E) = 1, a BCS-type DOS for the superconductor 
&s{E) = \E\/{\/E‘^ — A^) for \E\ > A and S>s{E) = 0 for |£'| < A, and a Breit-Wigner transmission function for the 
QD of the form[2] T{E) = (rir 2 )/(A£’^ + (Fi + r2)^/4) with /S.E = E — /iQD ± n • F^add (^ ^ N), which also accounts 
for gating of the single QD level by //qd = /^qd “ Vsd — electron filling with the addition energy 

F^add* 



FIG. S5: (a) Relative alignment of the electrochemical potentials of the superconductor /is, normal contact /xn and quantum 
dot /iQD for the bias applied to S. (b) G vs Fsd and Vg obtained from Eq. (1) with Fi = 10 /xeV, r 2 = 100 /xeV, Fadd = 7meV 
and A = 1 meV at T = 100 mK. 
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BOSON-ASSISTED QUASI-PARTICLE AND ANDREEV TUNNELING IN AN N-QD-S SYSTEM 


The subgap resonances parallel to the Coulomb diamond edges observed in our experiments (features IE and R in 
Fig. le of the main text) are due to resonant (R) Andreev tunneling (AT), with replicas due to boson-assisted, inelastic 
(IE) AT. Here we discuss in addition inelastic quasi-particle (qp) tunneling. [3] In this process a quasi-particle can 
be transferred from one lead to the other by the absorption or emission of a boson. To compare boson-assisted qp 
and Andreev tunneling, we briefly sketch the expected signatures of these processes in a charge stability diagram, 
see Fig. S6 for boson-assisted qp tunneling and Fig. S7 for resonant and inelastic Andreev tunneling. The energy 
diagrams to the right describe the processes at some characteristic gate voltages. 

Inelastic AT results in replicas parallel to the resonant Andreev tunneling line (R) with a spacing given by fico’b/2, 
while in boson-assisted qp tunneling the Coulomb blockade diamond edges are replicated with a spacing determined 
by A crucial difference between the two transport mechanisms is the arrangement of the boson-absorbing (red 
lines in diagrams) and boson-emitting resonances (blue lines). For boson-assisted qp tunneling shown in Fig. S6, the 
absorption replicas are inside the Coulomb blockaded regions and the emitting lines are above the Coulomb diamonds, 
similar to the structures observed for a QD contacted by normal metal contacts. [4] In contrast, the gate dependence 
of the Andreev processes is shown in Fig. S7. Resonant AT results in a resonance crossing the entire energy gap 
between the Coulomb diamond edges with a positive slope (green line). The boson-assisted replicas run parallel to 
this line, but in contrast to qp tunneling, the boson-emission and absorption lines merge where they meet the line 
connecting the Coulomb diamond edges with negative slope (see also Ref. 5). 

The nature of the observed resonances can now be determined based on the temperature dependence. With 
increasing temperature, the amplitude of an absorption line should increase due to the increased boson population, 
while the one of an emission line should decrease. In our experiments (see Fig. 3 main text and Fig. S13), we 
observe a decrease of conductance for the IE lines with increasing temperature and the appearance of absorption 
lines connecting to the emission lines, consistent with boson-assisted AT in the schematic of Fig. S7 and calculated 
in Ref. 5. We find no structures in the experiments related to the ones described in Fig. S6. Furthermore, for qp 
tunneling one could also expect broader subgap resonances and similar processes at the N-QD interface with different 
rates, generating another set of replicas sketched in Fig. S6 as weaker lines. 
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FIG. S6: Expected gate and bias dependence for boson-assisted quasi-particle tunneling in an N-QD-S system considering only 
single-boson processes. Blue lines denote emission, red lines absorption processes. Replicas of the gap edge can be either due to 
boson processes on the S side (positive slope) or on the N side (more transparent lines, negative slope). The energy diagrams 
to the right depict schematically the processes at specific gate and bias voltages. 
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FIG. S7: Expected gate and bias dependence of resonant (green line) and boson-assisted Andreev tunneling. Replicas of 
resonant Andreev tunneling can occur due to boson emission (blue lines) and boson absorption (red lines). 


EXPERIMENTAL SETUP AND SHIELDING 

The experiments were performed in an Oxford Kelvinox dilution refrigerator with an electron temperature of 
T ^ 100 — 150 mK. To minimize thermal radiation from warmer stages of the setup and the effects due to quasi¬ 
particles, [6] we filter the electrical lines using pi-filters at room temperature and a tape-worm filter between the IK 
pot and the mixing chamber. On the coldfinger directly above the sample holder additional pi-filters are employed 
to dampen also the IK radiation. The sample is shielded from the helium bath by a Faraday box and two additional 
copper shields. Since the electronic temperature is considerably lower than the critical temperature of the Nb strip 
(see section ) and the observed energy gap (see main text), the generation of quasi-particles is strongly suppressed. 
Furthermore, three independent cool-downs of the device in two different measurement setups allow us to exclude any 
spurious effects from the electronic setup (not shown). 


SAMPLE CHARACTERISTICS 
Characterization of QD 1 

Using individual charge stability diagrams as in Fig. Id of the main text and Fig. S8, we extract a leverarm of the 
global backgate of a ~ 90meV/V and an addition energy of F^add ~ 6.6 meV. The shift of the CB diamonds due to 
the superconductor dicussed in section and marked in Fig. S8 allows us to determine the superconducting transport 
gap at 5 = 0, Ao ~ 1.2 meV. In Fig. S8, we find conductance resonances running in parallel to the edge of the 
diamonds for |eUsD| > Aq. These are due to (electronic) excited states of the QD, which exhibit a Zeeman splitting in 
an external magnetic field [7] (see section ). From their roughly equidistant spacing in energy (marked in the figure by 
orange horizontal lines), we find a level spacing of dE « 1.2 meV, much larger than the spacing of the boson-replicas, 
5E ^ ficdb- The characteristic energies in our system can thus be summarized as /cT ^ F < ficdb SE ^ Aq F^add- 
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FIG. S8: Differential conductance G of QDl as function of the bias I/sd and the global backgate voltage Vbg at T = 110 mK 
and B = OT. The orange lines are used to read out the level spacing SE from the visible excited states. The white arrow 
indicates the visible superconducting transport gap Aq. Black dashed lines mark the shifted CB diamonds. 


No Cooper pair splitting or QD hybridization 


To exclude any effects due to the second QD of the device we measured G vs. sidegate and backgate voltage with 
the two QDs in series, see Fig. S9a. The resulting double-QD charge stability diagram in Fig. S9b confirms that the 
two QDs are well decoupled by the superconducting contact in between. In particular, we do not observe the typical 
honeycomb pattern or anti-crossings, but solely an increased conductance at the resonance crossing points. This 
suggests that the inter-dot coupling is much smaller than the individual QD life times, F 12 ^ rQDi,rQD 2 5 and that 
also capacitive cross-talk is negligible. In experiments with the two QDs in parallel, we do not find any conductance 
features that depend on both QDs, excluding effects due to Cooper pair splitting. 




-2.8 (V) -2.7 


FIG. S9: (a) Experimental setup to measure the two QDs in series. The superconductor is floating (not connected to ground), 
(b) Serial double-QD conductance map of sidegate 2 voltage vs backgate voltage at Vsd = OmV, B = 5T and T = 110mK. 


Additional analysis of the subgap features 

Figure SIO demonstrates the analysis of the subgap features and inelastic replicas (IE) of the resonant Andreev line 
(R) at 5 = OT and T = 110mK. Figure SlOa explicitly shows how we determine the position in energy of the IE 
replica lines relative to the resonant line R from conductance maps. The vertical cross-section through the middle of 
the features is used for the analysis of the conductance maxima Gmax and the full width half maximum (FWHM) of the 
IE replica lines. This cross-section also demonstrates the pronounced negative differential conductance and Fano-like 
linehape of the IE features, which agrees with the predictions in Ref. 5. From the analysis of the position and spacing 
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FIG. SIO: (a) Read-out of the energy Cn (and spacing) of the IE subgap features in a G vs backgate and bias voltage map. 
The crossing point between the negative slope [3- of the CB diamond and a black line along the conductance maximum of IE 
determines the position in energy of the feature relative to the resonant Andreev line (R). The IE4 replica (n = 4) at positive 
bias is pointed out as an example. The brown curve shows exemplary a vertical cross-section at the backgate voltage indicated 
by a dashed line, (b) Averaged position Cn vs. the IE replica number n, data averaged from 2 datasets for positive and negative 
bias. The red dashed line shows a linear fit Cn =< Se > -n, giving an averaged IE replica spacing of < Se >= 143fieV. (c) 
Spacing Se = Cn — en-i as function of IE replica state number n, data averaged over 2 datasets for positive and negative bias. 
The red line and area mark the averaged spacing < Se >= 143//eV and standard deviation determined in (c). (d) EWHM 
of the IE features as a function of n, for positive (blue circles) and negative bias (red squares). The line marks the value for 
r = 106//eV determined in the main text. 


in Fig. S10b,c we obtain an averaged spacing (over all n) between the IE lines of < 6e >= fico’b/2 = (143 ± 28) /reV, 
corresponding to an average boson energy of huj\^ = 286 fieV. 

Even if the IE replica energies scale roughly linear with n, some replicas significantly deviate from the average 
spacing (Fig. SlOc), possibly due to a different origin of the resonances, as discussed in the main text. Here, one also 
has to bear in mind that the position of Gmax might not completely agree with the position of the current maximum 
predicted for boson-assisted tunneling, [5] maybe partially explaining the deviations due to the overlap and slightly 
different widths of neighboring IE replica lines. To obtain an estimate for the width of the IE replica lines, we neglect 
the influence of negative differential conductance and determine the FWHM from Gmax to G = 0, plotted in Fig. 
SlOd as function of the IE number. We note that the FWHM stays roughly constant and smaller than the determined 
F from the fits, i.e. also FWHM < \/2r as predicted by theory. [5] We obtain an average FWHM of ^ 70/ieV for 
the IE lines. In summary, we find average values for the boson frequency of /b ~ 69 GHz, a FWHM ^ 17 GHz and 
a resonance Q factor of Q ^ /b/FWHM ^ 4. Similar to Ref. 4, one can study the maximum conductance of the IE 
replicas Gmax as a function of the replica line number (number of emitted bosons). We observe a varying Gmax (not 
shown), with the maximum for the negative bias e.g. at the n = 2 and n = 3 lines (visible in Fig. SlOa), reminiscent 
of a Franck-Condon blockade scenario. [4] Figure Sll demonstrates that neither the position nor the spacing of the IE 
replicas depends on the electronic charge state of the QD. The relative intensities of the replica lines are also similar, 
even if the total conductance of the replicas change for different CB resonances. 
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FIG. Sll: G as function of Vbg and Vsn over 3 Coulomb blockade resonances. The horizontal white lines mark the position 
Cn of the IE replica states read out in Fig. SIO, the line with slope /3_ is drawn for better orientation. The positions of the IE 
lines are independent of the charge state of the QD within experimental error. 


LINE SHAPE FITS AT FINITE BIAS 

To confirm that our fits in Fig. 2b of the main text at I^sd = 0 are independent of an overlapping Andreev 
bound state (ABS), we fit the expressions for a Breit-Wigner, the resonant AT and of a thermally broadened CB 
resonance to the conductance data at a small bias |Vsd| ^ Ao/e. Figure S12 shows such fits for a cross-section at 
VsD = —0.25 mV. The resonant Andreev tunneling line shape (red line) agrees very well with the data. From this fit, 
we extract Ti ^ 6.9//eV, F 2 ~ 68.3/ieV and F = Fi+F 2 ~ 75.2/xeV, in very good agreement with the values obtained 
in the main text for VsD = OmV. Due to the (weak) observed Andreev bound state (see main text), we tentatively 
attribute the larger tunnel coupling to the superconducting contact, i.e. Fs = r 2 . To control our fits in the normal 
state at larger magnetic fields, we also performed fits at 5 = OT and finite bias |Vsd| > Aq/c of the quasiparticle 
tunneling lines of the CB diamonds (not shown). Again the Breit-Wigner line shape agrees best with the data and 
we obtain very similar F-values compared to the fit at 5 = 5 T shown in the main text. 



FIG. S12: Fits of the resonant Andreev line R similar to the main text at B = OT and T = 110mK. The data (blue points) 
is a backgate sweep at Vsd = —0.25 mV. The data points to the right of the resonance are ignored (not shown) because of 
the inelastic AT features. The blue, red and dashed green lines are best hts according to the Breit-Wigner, resonant AT and 
thermally broadenend GB line shape, respectively (see main text). 


EVOLUTION WITH TEMPERATURE 

Figure S13a shows G as a function of Vsd and Vbg foi* different temperatures T = 110 mK, T = 500 mK and 
T = 950 mK. At higher temperatures, kT ^ /kJb/2, the conductance maximum of the resonant Andreev line and the 
emitted-boson replicas start decreasing and boson-absorbed sidebands appear. This characteristics is already evident 
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FIG. S13: (a) G vs Vbg and Vqb at T = 110 mK, T = 500 mK and T = 950 mK. The dark dashed lines mark the position of the 
studied horizontal cross-sections at hxed I^sd, shown exemplary in the T = 110 mK map. The shaded area at Ksd = —0.5 mV 
visualizes the data used for averaging. White arrows in the T = 950 mK map point at boson-sidebands due to boson-absorption, 
(b) Similar conductance map as in (a) at T 5K. Even at this high temperature similar lines and subgap features are still 
visible, marked by white arrows, (c) Extracted EWHM vs temperature for the n — 1...4 IE replicas for positive and negative 
bias. The EWHM is constant for this temperature range and slightly smaller than T = 106 /xeV extracted for the QD life time 
broadening (from the main text, dark blue line). 


from the raw data (white arrows in Fig. S13a). We study this behavior also in cross-sections at fixed Vsd values, 
indicated in Fig. S13a. To average out noise occuring at higher temperatures, 12 cross-sections around the desired 
Vsd value were averaged, including the compensation of the horizontal shifts due to gating. Other cross-sections (not 
shown) show the same features as the Vsd = —0.5 mV cross-section discussed in the main text. From a similar analysis 
as in section and Fig. SIO we can conclude that both the position and the spacing of the IE replicas stays constant 
within the studied temperature range. Fig. S13c shows the extracted FWHM of the IE replica lines as a function of 
temperature. The FWHM of the resonant line R (not shown) and the replicas stay roughly constant and FWHM < F, 
in agreement with the prediction [5] FWHM < \/2F independent of T, further supporting our interpretation. Fig. 
S13b shows a conductance map recorded at T ^ 5K < Tqi = 7.7K (cf. section ). Even if most features and the 
superconducting gap are already smeared with temperature, one still observes pronounced lines and features in the 
subgap region running parallel to the positive slope of the CB diamonds, marked by white arrows in the figure. 


EVOLUTION WITH MAGNETIC FIELD 


Figure S14 shows G vs. Vbg and Vsd for a series of external magnetic fields applied perpendicular to the sample 
plane. From such conductance maps we extract for each field in the main text the following quantities: (i) the 
transport gap A from the tips of the shifted diamonds as discussed in section , (ii) the Zeeman splitting of the ground 
and orbital excited states of the QD, indicated in Fig. S14d, (iii) the 5-field dependence of resonant Andreev tunneling 
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FIG. S14: Example data for the B-field analysis: G vs Ebg and Vsd colormaps for (a) B = 0.25 T, (b) B = 0.5 T, (c) B = 1.5 T 
and (d) B = 5T. The magnetic field is applied perpendicular to the sample plane and CNT. The vertical dashed line in (a) 
indicates the position of the measurement in Fig. 4a of the main text. In (c) and (d), the marked and AE^^^ correspond 

to the Zeeman splitting of the ground and orbital excited states for positive and negative bias, respectively. 


from horizontal cross-sections at Vsd = 0 mV, (iv) the position, spacing, conductance maximum Gmax and FWHM of 
IE replica lines from an analysis similar to section . We now discuss some additional results in more detail. 

Figure S15a and b show the analysis of the Zeeman splitting AEz = of the ground and first orbital excited 

state of QDl. We obtain a g-factor of g ~ 1.9. Figure S16 shows the magnetic field dependence of the resonant 
Andreev tunneling peak R, for (a) the conductance maximum Gmax and (b) the extracted Fi from individual fits to 




FIG. S15: Zeeman splitting AE^^ (a) and AE^^^ (b) of the ground and first orbital excited state for positive and negative 
bias, respectively, read-out from the B-held colormaps as demonstrated in Fig. S14. We obtain a g-factor of roughly g ^ 1.9 
from the fits AEz — ggsB. 
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(a) 




FIG. S16: B dependence of resonant AT. We study horizontal cross-sections at Vsd = OmV as function of extracted 
from individual colormaps. (a) conductance maximum Gmax of the resonant line (R) vs B. (b) Fi extracted from fits to the 
expressions for resonant AT (red circles) and the Breit-Wigner line shape (blue squares). The shaded red area indicates the Fi 
values for fits in good agreement with the data. 


the resonant AT and the Breit-Wigner (BW) line shape (see main text). Fi is the most sensitive parameter in these 
fits. In both plots, we observe a clear transition (marked on top of the plots with a colorscale) with contributions 
from resonant AT only (red, up to B ^ IT) to an intermediate region where both processes coexist {B ^ 1 — 3T) 
to contributions from normal electron tunneling alone (blue, B « 3.5 — 5T). This is immediately evident when one 
looks at the range of Fi for acceptable deviations from the high-held BW values in Fig. S16b. The values for small 
B extracted from resonant AT and for large B from normal tunneling agree very well. Similarly, one observes an 
increase in the conductance from Gmax ^ 0.06 x 2e^/h to Gmax ^ 0.22 x 2e^/h in the intermediate region where 
both processes coexist, followed by a decrease of conductance to Gmax ^ 0.16 x 2e^/h at high helds. From the 
conductance maps, we can also extract the visible transport gap A as function of 5, shown in Fig. S17a. We hnd 
that A is suppressed at 5 ^ 1.5 T. This held corresponds to the transition region from resonant AT to BW. The 
detailed characteristics of A seems at hrst somewhat surprising. A is reduced fast up to 5 ^ 0.3 T, then stays roughly 
constant up to 5 ^ 1.2 T, before it monotonuously tends to zero ai B ^ 1.5 T. In the region between 0.3 and 1.2 T, 
we observe non-monotonuous, systematic huctuations around a constant value. The positions of this opening-closing 
behavior of A are quite reproducible in several repeated measurements and also evident in the Vsd vs B map shown 
in Fig. 4a of the main text. In addition we observe a hysteresis in A between up- and downsweeps of B (not shown). 




FIG. S17: (a) Transport gap A vs extracted from individual colormaps similar to Fig. S14. The red dashed line is a guide 
to the eye. (b) Extracted FWHM vs B for the n = 1...4 IE replicas for positive and negative bias. The EWHM stays constant 
as function of B and smaller than F = 106/xeV (from the main text, dark blue line). 
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We ascribe this behavior to vortices in the superconducting type (II) phase of the narrow Nb strip. Indeed, a first 
critical field Bqi ^ 0.3 T, marking the initial decay of A and the onset of the hysteresis and ‘jumps’, would be in 
reasonable agreement with literature values for thin Nb films. [8] For larger magnetic fields, vortices could enter the 
region of the narrow Nb strip close to the CNT. We speculate that, depending on the position and distance of these 
vortices relative to the ‘sensing’ CNT QD region the visible gap A could be larger or smaller, thus explaining the 
sudden jumps in A with a vortex rearrangement. As an extension to the analysis of the position, amplitude and 
spacing of the replica lines in a magnetic field in the main text. Fig. SI 7b shows the extracted FWHM of the IE 
replica lines as function of B. The FWHM stays constant as function of B and FWHM < F holds. 


A SECOND SAMPLE 


Figure SIS shows supporting data with very similar findings for QD2 of the same device and for a second sample 
fabricated in the same geometry. In both cases, one clearly sees straight lines running parallel to the edge of the CB 
diamond as discussed in detail for the first sample. 


(a) G (2e2/h) 



-2.38 -2.36 



FIG. SIS: C vs the backgate voltage and bias Fsd for (a) QD2 of the same device (b) an independent sample with the same 
geometry, recorded at T = 60mK. In (a), weak conductance lines running parallel to the CB diamond are labeled with white 
arrows. In (b), an Andreev bound state (ABS, orange arrows) is more pronounced, but a set of lines running parallel to one 
edge of the diamond are still visible. The resonant Andreev line R and inelastic replicas IE can be clearly identified here. 


CHARACTERIZATION OF THE SUPERCONDUCTING NIOBIUM STRIP 

Figure SI9 shows the resistance Rstrip of the Nb strip acting as S-contact in our devices as function of (a) the 
external magnetic field, and (b) the temperature, measured in a two terminal geometry with a cryostat line resistance 
of (a) 415 U and (b) 2.5 U. The magnetic field dependence in (a) shows a broad transition to the normal state at 
^ 3 — 4.5 T, which coincides with the transition to the pure Breit-Wigner line shape discussed in the main text (cf. 
section ). In (b) we find two distinct superconducting transition temperatures, most likely due to different widths in 
the design of the superconductor strip. We ascribe the lower Tqi ^ 7.7 K to the narrow part of the Nb strip in direct 
contact to the CNT. In the experiments presented in the main text the sample temperature is well below Tci, so that 
we can neglect a temperature dependence of the superconducting energy gap. 
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FIG. S19: Resistance of the Nb strip measured in a two terminal configuration as a function of (a) the magnetic held at 
T = 110 mK base temperature, and (b) the temperature in a variable temperature insert (VTI). The two arrows mark two 
distinct superconducting transition temperatures Tci ~ 7.7 K and Tc 2 ~ 8.7 K. 
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